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ABSTRACT 


The underlying purpose of this research was to understand 
and predict the response of upper ocean boundary layer 
mixing to high frequency oscillations in the atmospheric 
Mmereing attributable to the diurnal heat flux cycle, 
unsteadiness in the wind, and other short term changes. To 
accomplish this task, a non-stationary, one-dimensional 
bulk model of the mixed layer originally proposed by 
Garwood (1976, 1977) is generalized by deriving a new 
equation for the entrainment buoyancy flux which includes 
the unsteady term. To examine the importance of the un- 
steady term, a reduced form of the turbulent kinetic energy 
equation is solved numerically. The results predict the 
high frequency cutoff above which the unsteady term should 
not be neglected. The quasi-steady state assumption is 
feemrace Only for low frequency forcing up to this cutoff. 
There was also some dependence upon the amplitude of the 
imposed forcing. To investigate the asymptotic properties 
of the mixed layer model, the full form of the Garwood 
model is solved numerically and compared with the solutions 
to the de Szoeke and Rhines model. The results of this 
research provide a theoretical basis for realistically 
applying the mixed layer models to ocean boundary layer 


problems on all time scales. 
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I. INTRODUCTION 


Peer URPOSE OF THE STUDY 

The objective of this research was to study the one- 
dimensional mixed layer model of Garwood (1976, 1977) in 
terms of a possible generalization of the extrainment 
buoyancy flux equation and to examine the responses of the 
mixed layer turbulent energy intensity to the different 
amplitude and frequencies for wind-shear turbulence prod- 
Potion. Because of the frequency dependence, the steady 


a<E>_ 


ete ( —=-~=0 , where <E> is the turbulent kinetic energy 


averaged over the mixed layer depth) or quasi-steady state 


a <E> 
at 


but it is not possible in other ranges. Using a constant 


az 0 ) assumption is good in some frequency ranges 


dissipation time scale, Garwood (1976) suggested that if 
the quasi-steady state assumption is made to facilitate the 
Solution without filtering the surface boundary conditions, 
an incorrect high frequency response will not only be pres- 
ent but may bias the mean mixing trend. 

erurEner ObJective of this study is to critically 
examine the asymptotic time regimes of the Garwood model 
in comparison with those of the de Szgoeke and Rhines (1976) 
model. These two tasks will provide a theoretical basis for 
realistically applying the mixed layer model to ocean 


boundary layer problems on all time scales. 





Be METHOD 

The equations of the Garwood (1976, 1977) model are 
listed and their derivation is discussed briefly. The 
de Szoeke and Rhines (1976) work is also reviewed in 
order to provide the background information necessary to 
compare the asymptotic properties of the two different 
models. To generalize the Garwood entrainment buoyancy 
flux equation, the unsteady term of the turbulent kinetic 
energy equation at the base of the mixed layer, which was 
usually considered to be negligible oy earlier investigators, 
is parameterized. For the study of frequency-dependent 
response to atmospheric forcing, the method applied by 
Garwood (1976) in which he used a reduced and nondimensional- 
ized form of the turbulent kinetic energy equation, and 
solved analytically by using a constant dissipation time 
Scale, is reviewed. This same form of the equation is 
Solved numerically by using a non-constant dissipation time 
scale, h<E»* , provided by Garwood (1976), where kh is the 
mexead layer depth. In this case, the solutions for non- 
dimensionalized frequencies Wh/Us=20071, 2nN , .277 , 
and .Qd2%T , where W is the angular frequency of the 
Memmalized wind forcing and u, is the maximum friction 
memecity, and amplitudes of fluctuations a= 0.1. 0.5 and 
1 for the wind production are compared. 

To study the asymptotic properties, the Garwood model 


is first solved numerically for different stabilities and 
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its results are compared with those of de Szoeke and Rhines 
by using equivalent calibration constants, boundary con- 


mepr1ons and initial conditions. 


C. HISTORICAL LITERATURE REVIEW 

The vertical fluxes of heat, salt and momentum across 
the sea-air interface are the source of almost all oceanic 
motions. In the fully turbulent oceanic mixed layer 
bounded by the sea-air interface above and by the dynamically 
stable water mass below, the vertical fluxes are large. 
Below the mixed layer they are usually negligibly small so 
that one can decouple the mixed layer from the underlying 
Stable, quiescent water mass. The immediate local reaction 
of this mixed layer to those fluxes result in a homogeneous 
feeser column, 1.e., vertical uniformity in the mean velocity 
and density. This adjustment of the density and velocity 
structure of the surface layers of the ocean to variable 
fluxes has been the subject of a large number of studies 
Since Ekman (1905)'s treatise where he originated the 
@eneept of a depth, d = 7. 5isec) Mo, (W is the wind velocity 
and 9 is latitude), of frictional resistance for the upper 
Section of a wind stressed ocean and suggested that any 
Surface mixed layer by the action of wind has the same 
order of this depth. This depth comes from the mathematical 
solution to the steady state horizontal momentum equation 
in which Reynold's stress is related to the mean shear by 


a constant eddy viscosity, K . Rossby and Montgomery (1935) 


ILS 





Bemmbed OuL that the depth of a surface drift current layer, 
het,/f (fis Coriolis parameter), and Ekman depth are not 
necessarily comparable. Both depths are derived by con- 
sidering only the momentum budget, neglecting the effects 

of buoyancy and the mechanical energy budget. 

Munk and Anderson (1948) feralowed a simultaneous 
solution of steady-state heat and momentum profiles by 
using a viscosity and an eddy conductivity that depend upon 
the gradient Richardson number. They did not recognize the 
presence of a sharp interface at the bottom of the mixed 
layer and therefore their result is closer to Ekman's 
Solution than to the physical reality. Kitaigorodsky (1960) 
concluded by dimensional analysis that the mixed layer depth 
tec De proportional to the Obukhov length scale, L. The 
basic flaw in his steady state model is the assumption that 
the entire ocean mixed layer is analogous to the constant- 
flux atmospheric surface layer, because with his constant 
flux at the surface the mixed layer temperature and depth 
cannot both remain unchanged. 

Kraus and Rooth (1961) developed a plausible steady state 
model based primarily upon the buoyancy equation. In their 
model, steady state is achieved by balancing solar radiation 
with a net surface heat loss due to evaporation, conduction 
and long wave back radiation. The unstable density profile 
above the compensation depth is a convective source of the 
turbulent kinetic energy for mixing. Recognizing the limi- 


tation in Kraus and Rooth model--no provision for a possible 
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downward surface heat flux, no account of mechanical 
@eoadction Of turbulent kinetic energy and the steady state 
constraint--Kraus and Turner (1967) improved the one- 
dimensional model by considering the turbulent kinetic 
energy budget utilizing the two separate equations, the 
heat equation and the mechanical energy equation. Because 
the frictional generation of heat is negligible, the 
vertically integrated heat equation provides a relationship 
for the conservation of potential energy. They parameter- 
im7ea the mechanical production rate in terms of the friction 
velocity but neglected the viscous rate of dissipation and 
the effect of entrainment shear production. From the 
examination of observations, Turner (1969) deduced that a 
substantial fraction of the part of the work done by the 
wind which goes into the drift current is eventually used 
to deepen the surface layer. 

Miropol'skiy (1970) and Denman (1973) assumed that 
Mfesipavtion is a fixed fraction of mechanical production 
and the remaining turbulent kinetic energy goes to downward 
buoyancy flux. Pollard, Rhines and Thompson (1973) apply 
the slab approach to the oceanic mixed layer but they com- 
plete the entrainment problem with a different mechanical 
energy requirement. In this model, the energy for the 
entrainment is derived directly from the mean flow and the 
intensity of the turbulence is not considered to have an 
active role in the entrainment. Niiler (1974) in his three 


layer model which is a combination of Kraus and Turner, and 


AY 





Pollard et al considered that both turbulent kinetic 
energy and the mean kinetic energy are important in the 
entrainment neeen: * Although he added the entrainment 
production term, the surface production minus dissipation 
was still parameterized in terms of ux . Therefore, 
dissipation was not permitted to adjust to other varying 
atmospheric conditions or to entrainment shear production. 
A more realistic parameterization for dissipation is not 
possible without actually computing the mixed layer turbu- 
lent kinetic energy intensity <B>. It cannot be done with 


only the surface scale ux. 


fee REVIEW OF THE RECENT LITERATURE 

1. TDhe Garwood Model 

Garwood (1976, 1977) developed an ocean mixed layer 

model using the Navier-Stokes equation of motion with the 
geostrophic component eliminated, the continuity equation 
in incompressible water, the heat equation from the first 
law of thermodynamics, the conservation of salt equation, 
and a linearized equation of state. A summary of the 
modeled equations follows: 


entrainment buoyancy flux, 


—~ Ew (=h) - MF <5? ie) 
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budget for horizontal components of turbulent 
kinetic energy, 


( d 12 2 L bwi(-h) _ EV. ’ =v 
x ie (hsurey’ CBN) S401 Te eee «OE aN >< =? 
ee E> ial 


budget for vertical component of turbulent kinetic 
energy, 


4_ (hh cam>)s HEWEH 


reLAPey, = = 
> hbw'(o) + (<B> ~30W™>) CE 
dt 2 


nN |— 


- £(¢ 574 $h)<E> (1-3) 


conservation of mean buoyancy and mean momentum 


n.d<b> wabw'(o) 4 OWi(-h) 4+ ee J, Qdz 





(1-4) 
db (h¢ey)z - Tw Co) =F F<2>h Gie5) 
\Umpecondi tions gat bottom of mixed layer 
-Fw’ (-h)= ab SEA (1-6) 
OW (-h)= ot FE (1-7) 

fee the following notation: 

Uy = Uy + U: i Ea ne re ) dx dy 
Gc zu +iV Q wear ie 

Es Wt + v'74+ wi? fis weer etaes 

ug 2 | 7% ee a. 
Peenie.-€) 2/ @, 

4243 

ie net 
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To derive the equation (1-1), the equation of the 
turbulent kinetic energy per unit mass was derived from 


the equation of motion using the Boussinesque approximation, 


ae: 2 (er a (pe 77 Ou; __ 8 ieee 
Z(E)+ 5 OX; (Su;) + Pa, (ui; E) 4 U; UY; Be See oe, 

P’ ou} ou; ae y Bui ou; aS ; = 
te 3x, ax; ae a or ui SS) » 3x DX; + Ub 453 (1-8) 


Using the boundary layer approximation and the continuity 
equation, and assuming that viscous diffusion and dissi- 
pation of the mean kinetic energy are negligible, the 


equation (1-8) was written: 


ec ee ee + Vw 20] GW -¢ 


where 


coy 
ake (uy auf) 

With the assumption that the turbulence of the overlying 

mixed layer provides the energy needed to destabilize and 

Seode the underlying stable water mass. The local (at z=-h) 

turbulent kinetic energy budget is the basis of the entrain- 

ment hypothesis. Assuming a quasi-steady state of z=-h, 

and that within an active entrainment zone, the most 

Significant source of energy for mixing is the convergence 

of flux of turbulent kinetic energy —- 2 [ws +E )| , equation 


(1-9) was reduced to: 


Bwh)= 2[w(E+E)]L, 


Craik.) 


20 





Estimating the time scale, tf. required to transport some 
of the turbulent kinetic energy «E>, to the viscinity of the 


entraining interface, 


_ = i( E ee _ “E> ol) 
Sz lw ( e+e), Te 


and assuming that the mixed layer depth is proportional to 
distance over which turbulent kinetic energy must be 
transported by the vertical component of turbulent velocity, 
iw?sx?, G was taken to pehKw and therefore from (1-10) 
and (1-11). 

Lw'?> <=> 


bw’ (-h) ~ 4 Geb 


Memeec the equation (1-2) and (1-3), the equation (1-9) 
Meoeitirst divided into the horizontal and vertical components 


by assuming an isotropic dissipative structtrre and recalling 











that Ese = 0 is included in the equation (1-9), 

1 otutsvi. 3 w EE) | aw 24 yaw Pou , pr av’ _2 eke 
— oth 2 2 eal ad az + & ox 7 © ay ac 
and 

LW ea lw (Wek) 4 Pow atw- te Geils) 
z ot az 2 e, C, az 3 a 


the equation (1-12) was integrated over the range from z=-h-& 


Mem2=0 (Figure 1) 


ees a dz =- i = [ w’ (eee) ] dz- ie [aw 24 aviw! at wy Jy 
4(° £ 1G (34 + Tg je 2 ow 7 ets 
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Applying Leibnitz's rule and using the assumption (U“7+V") (-h-iE0, 


emo (uv) ss : —_—_——> = 
‘G >t dz = df (ua $v?) dz d(h&u +v'*)) (1-15) 


Defining a dissipation time scale as a Edz = 8 Mince 
the assumption of locally isotropic turbulence and con- 
Sidering that Ge = (Ga + om in which Te, is the 
convective time scale of large eddies (proportional tohkE™ 


and @2z 1s the time scale of planetary rotation 


(proportion to fa ), therefore te = ER" 4 f and 


7 . é dz 


The vertical integration of the pressure redistribution term 


Ae iingle Evie (1-16) 


fmeman important source or sink term for the eave 
turbulent kinetic energy budget even though ra are 
eC 3x; 


Following Rotta (1951), Lumley and Khajeh Nouri ee 


equivalent bulk formulation was 





9 you; St dz = cEY” (> S51e URS) (| alte), 
4-5 & 3X; 


6 ‘ 7 ! ——— 
ae eo ~- c EY* (E> - 34° i 
Therefore, - 6, ( Sy gad av) dz 2 < E> (> 3¢ >) (1-18) 


Following Kraus and Turner (1967), Denman (1973) and Niiler 


(1975), the shear production term was parameterized: 


¢ “Ta? al —_ = Cc |ocl dh - 
S. .(aw 1S + Wwe) dz sue OL + - (219) 
where|fG| was the "excess" surface mean velocity in the 


direction of the wind stress. In this case, the inhomogeneity 


of the mean velocity field cannot be neglected. 


Ze 





é 
-(° (2 w(ByE))] dzem'u 
~~§ oz e, 2 
Tft§&l is proportional tou,, then (1-18) may be combined 
with the parameterized net input from breaking waves, less 


the loss to radiating internal waves; 


4 


“52 [Beto oB0 0m Ba 7 BE] derma BB (2-20 
“h-6 e 


Peeome (1-15), (1-16), (1-18) and (1-20), (1-14) was written as: 


L har evityys m Ug SEL GA ~EN4 E> -3CwD) - 2 CES*4 Fh) ED (1-2) 


Similarly integrating (1-13) over the range from z=-h-é to Z=o, 


o d _(°? 74 BP ow’ - w(“e+ek 
Z), gins wr) dz = J) CBW + & SF ) dz War +&)}, 


WEL 5 E dz 


ion 
and assuming Wo at the surface and wel f negligible, 
the vertical component of turbulent kinetic energy budget 


was parameterized by using (1-16) and (1-17); 


(be wr) 2 dbp Bw (0) + Bw) ] Ey (Ey -3 60>) (1-3) 
~L(cey*4 fh) <ED 


From the heat equation, conservation of salt equation 


aL dW 
2 at 


and equation of state, the buoyancy equation was formed. 
Averaging this buoyancy equation and using the boundary 


layer approximation, the mean buoyancy equation was obtained: 


2b = aw’ 1 <2 
sa + Be (es2ib) 
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To get the equation (1-4), (1-21) was first integrated 


over the range from Z=-h-9d to Z=0 


o a bw! J re) _ 
Fe: at 22 dz = “lt "or foze dz + 58 (Adz (1-22) 
Using the Leibnitz's rule and taking the limit 2-70 ; 
Gae term on the left hand side can be written; 
ab ee d(-h-5) 
Jans 3¢ eZ oan bdz -[-bCh-5)] “aE ott 


—— dh 
~ Se (h< b>) —6Ch-7) a 


Assuming bw’ (-h-5)=0 and using Abz¢b> -b(-h-$) , (1-21) was 
written; 


p dee _ a5 dh Bw (o) + 429" dz (1-2l1) 


To get the jump condition (1-6), the equation (1-21) was 
integrated over the range from Z=-h to 220 and from Z to 
Z20 , and these two integrated equations were solved 


Simultaneously to eliminate ach » Then the result was 


pw’ (z) 2 BW(o) (1+ B) -(Bwi Ch) + 3% 7. Qdz)= 


ae) Adz 


Smiiarly the integration of the equation (1-21) over the 


(1125) 


range from Z2-h-§ to z=-h with dropping the negligible solar 
radiation term and assuming Vw’ (-h-3)=o0 and approaches 


zero led to 


2S = “BW (-h) (1-26) 
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By applying the Leibnitz's rule, 


+ ak py dh) tps deh) ; 
op ot: ys ie [ Bch) Se - ~b¢h-) at (a2) 


and by assuming lim d(S5<bY) — 0 and ye dd =06 


wr 0 dt 
-h ; dh _ at dh 
jim Fg Be dz=[BCm-bCh-A)] G = ab RA (1-28) 


Therefore from (1-26) and (1-28), the jump condition (1-6) 


meas obtained 


-Gw' Ch) = ab (1-6) 


and using this relation the equation (1-24) can be written 
in the form of the equation (1-4). 

Using the boundary layer approximation and dropping 
the negligible viscous terms, the equations of mean 


motion were combined into a complex form by using @=UFti¥ 


ac mafz. sow 


> => (1 =29)) 


To get the equation (1-5), (1-29) was integrated vertically 


over the range from Z2=-h-d to Z=0 with the assumption of 

CW! (-h-¢) =0 ; 

2 (h¢e>) = ~ifheey - FW! (0) (1-5) 
The jump condition (1-7) was obtained in a similar way as 


(1-6). 
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The most important aspect of this model was the 
ability to predict the year-around evolution of the mixed 
layer depth together with bulk properties. This model 
provides not only structural prediction but also predicts 
a changing layer depth because the entrainment model pro- 
vides boundary conditions at the moving density interface. 
In addition to the non-linear dependency upon stability, 
he found a Rossby number - R= Ux/(fh) - dependence for the 
entrainment rate in his nondimensional solution. This 
makes possible a cyclical steady state for the boundary 
layer without requiring unrealistic values of upwelling or 
lateral advection Pia wong GeEM in besrabL7on. 

On short time scales of the order of a few days, 
the upper ocean, even at higher latitudes, exhibits 
Significant baroclinec activity. The importance of shorter- 
meried fluctuations in the surface buoyancy flux in 
modulating the long-term response has demonstrated the 
need to know the typical daily heating/cooling cycle, solar 
radiation, evaporation, conduction and back radiation as a 
function of season and geographical coordinates. Some of 
the most fruitful applications of this model lie in the 
link-up with other models, for example, joining with an 
ocean circulation model. 

In summary Garwood suggested in his model that first, 
the planetary rotation 1s assumed to influence the dissi- 
pation for deeper mixed layers and enables a cyclical 


Steady state on an annual basis. Secondly, the rate of 
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entrainment for the stable regime is not accurately 
reflected by a linear extrapolation of the unstable 
Situations. This is particularly important in modeling 
the ocean boundary layer. Unlike the atmospheric 
boundary case, most of the solar radiation does not 
penetrate the layer. Therefore, downward turbulent heat 
flux in the oceanic boundary layer is as important as the 
upward flux during the course of both diurnal and annual 
cycles. The non-linearity of interface Beer ecnent 
tendency parameter, which is greatest for stable surface 
boundary conditions, results in a modulation of the long- 
term trend of mixed layer depth by the diurnal component 
Semsuriace heat flux. Thirdly, in this model buoyant 
production is somewhat more efficient than shear production 
as a source of energy for vertical mixing because of its 
unique effect on the vertical component of the turbulent 


melocity. 
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2. The de Szoeke and Rhines Model 





The de Szoeke and Rhines (1976) study of initial 
mixed layer deepening is reviewed here because the results 
will be critically analyzed in a comparison with the 
Garwood model. The solution to impulsive initial deepening 
is most relevant to this research on high frequency response 
of the upper ocean. Following Niiler (1975), de Szoeke and 
Rhines (1976) investigated the turbulent energy balance of 
a wind-mixed layer (in the absence of solar heating) by 
adding a new term which represents the energy to spin-up 
the level of turbulence in the increment of the mixed 


layer dh in time dt. The equation is 


{ vet Zul (1~ cos ft) Fa (Ge 7 | dh = 2m, u2h (1-30) 


ct 


(A) (B) (Cc ) (D ) 


where term (A) represents the entrainment damping, term (B) 
represents the shear production at the base of mixed layer, 
term (C) represents the energy needed to spin-up the level of 
turbulence and term (D) is the turbulence production minus 
@tsesipation, in which it 1S assumed that the dissipation is 
pened fraction of turbulent production. 

The results (using constants and parameter m,=G =! , 
Up=i “m/sec , Laio0* sect , Pan2= 0% sec* ) were 
Mmemetor slightly later than the initial instant, since 
Cos $t = 1~ Ea , the shear production, term (B), can be 


reduced tou, tt? and 1s negligible compared with other terms. 
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Assuming a negligible entrainment damping at this early 
stage, the balance (C)~(D) holds for small t. Initially 

it is a linear deepening, i.e., dh cc Us and this lasts 
for a fraction of a Brunt-Vaisala period (~100 sec). It 
attains a depth of 2~7m while the level of turbulence in 
the incipient layer spins up. The term (C) here represents 
a Significant demand on the available turbulent energy only 
during the first Brunt-Vaisdla period after initiating the 
mixing. After that, it may be safely neglected. 

mee For large t (about 100 sec to 1 hr), a balance holds 
between (A) and (D). After it starts deepening by turbulent 
energy, the entrainment damping term grows. A +t” deepening 
emerges as the turbulent work done by the wind at the 
surface erodes heavy fluid from below and mixes it 

uniformly through the surface layer, raising the potential 
energy. This one-third power deepening proceeds from the 
Same mechanism that was suggested by Kraus and Turner 

(1967) 

jm) For still larger t (about 1 hr to 12 hr), a balance 
(A)~-(B) holds approximately The mean flow has been 
accelerating until after 20 to 120 mins ata depth of 5~—20m, 
the production of turbulent energy by the shear at the base 
of the mixed layer dominate other sources and a t”%* deepening 
takes over. This h=t¥’¥? was proposed by Pollard, Rhines and 
Thompson (1972). 

iv) For later times, the Coriolis force swings the mean 


motion vector away from the surface stress direction, 


ee 





reducing the entrainment shear production (term B) below 
the rate Meee assumed above. As t increases substantially 
beyond t=7/$ (corresponding to the lapse of one-half 
pendulum day). An approximate balance between (A) and (D) 


is established again, leading to tt’? development. 
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II. GENERALIZATIONS OF ENTRAINMENT BUOYANCY 
FLUX EQUATION OF THE GARWOOD MODEL 
Garwood (1976, 1977) developed the entrainment equation 

from the total turbulent kinetic energy budget by assuming 
that the turbulence of the overlying mixed layer provides 
the energy needed to destabilize and erode the underlying 
stable water mass. The total turbulent kinetic energy 
equation (1-9) at z=-h is: 

3E 2 ou Ov 
5 ial =. ee (fE P? = hime. + ow 
2 at '-h az LW (E+E )), [ Pw =p oz J-h 


+ Bw}, -e [4 e) 


The importance of the unsteady term in the atmospheric 
energy budget was demonstrated by Zilitinkevich (1975) in 
the case of weak stratification, and also suggested by 
Garwood (1976) for the oceanic planetary boundary layer. A 
more detailed study of the importance of this term will be 
in the next chapter. For the order of magnitude of the 
Shear production term, Garwood (1977) found that shear 
production is a fixed fraction of buoyant damping in the 
entrainment zone. Zeman and Tennekes (1977) also suggested 
in their atmospheric planetary boundary layer model that shear 
production is absent if the wind shear across the inversion 
base is zero. Even if there is substantial shear, the 
effect of this term appears likely to be relatively small 


(Mahrt and Lenschow, 1976), because mechanical production 
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of turbulence tends to be largely balanced locally by 
viscous dissipation through the energy cascade mechanism 
(Tennekes and Lumley, 1972). Since the position ( Z=-h _) 
of the entrainment zone of the oceanic planetary boundary 
layer is defined by the location where the turbulent 
kinetic energy vanishes, local dissipation is also assumed 
to be a fixed fraction of the flux convergence yeion, Saennille 
it be significant in the energy budget. Therefore, follow- 
ing Garwood (1977), a balance holds among the flux 
convergence term, buoyancy damping term and unsteady term. 


The equation (2-1) is reduced to: 


ioe!) . Wh Ase a = 
2 ln el (Ste) +t PL, si 


Since the boundary erodes and retreats depending on the 

level of the intensity of the turbulent kinetic energy, 

it may be assumed that the level of turbulent kinetic 

enerey at the moving interface remains unchanged in the first 


approximation. Therefore, 


dt leh at I-h az at eo © 


or (2-3) 


25 | ~_ 9E 22 
atl-h ~ OZ ati-h 


Assuming that the turbulent kinetic energy decreases 
linearly with depth in the mixed layer and using the finite 


difference approximation, (2-3) can be written: 
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E E a(-h €)-0 = 
ee soe So. (a8) Ot, 2 Bh (2-H) 
ot i-h 4z2 ot “3 -Gh) ot h ot 
where <E» is the average value through the mixed layer. 

; ; : h 
Employing the parameterization (1-11), Te = tany , and the 
jump condition (1-6), 
| 2¢— ah <B> on, Er QW) 4 oa, 
ae st ayy a aay ra bw’ (-h) 


and this forms 


r raw 
Cw! (~h) = E> Cwey” 2-5) 
h +(E2/AaBb 


This more general equation of entrainment buoyancy flux 
meplaces the equation (1-1) of the section I.D.1 for the 
final closure. 

Therefore, the entrainment buoyancy flux, 
and the rate of layer deepening from (1-6) are functions 
of the bulk values, or overall layer averages, of both 
the total turbulent kinetic energy and the relative dis- 
tribution of this energy between the horizontal and vertical 
components. The advantages are two-fold: 

i) The rate of deepening is an integral property--not 
just a function of the local properties of the turbulence 
field, which may be vanishingly small at Z=-h inspite of 


a Significant rate of deepening. 
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ii) The relative distribution of turbulent energy among 

the three geometric components is taken into account. 

Being able to include the root-mean-square vertical 
velocity, CW'2>%, is particularly important in realistically 
treating layer "retreat", =m <0. 

In the case where WaBb ED » equation (2-5) 
meeurns to (1-1) of the original model. On the other hand, 
in the case <E)> » hab which is the case of weak 
stratification comparing to turbulent intensity (2-5) 
reduces to a = (way? Meeirso tse it the stratification 
is close to the neutral lapse rate, the entrainment proceeds 
at a rate proportional to the root mean square vertical 
velocity, <w*"> . The comparable situation in the Zeman 
and Tennekes model is that §w’ tends to zero as the lapse 
rate, ~ , tends to 0 (adiabatic conditions). Thus the 
energy budget in the absence of shear term and dissipation 
terms becomes a balance between the unsteady term and 
fiux. convergence term i1.ée., “t Cc Ow where Ow is 


the average standard deviation of vertical velocity within 


the mixed layer and will scale with {wre : 
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III. DISSIPATION AND FREQUENCY RESPONSE 


A. DISSIPATION AND ITS TIME SCALE 

Turbulence is always dissipative and this dissipation 
occurs through viscous shear stress doing deformation work, 
which increases the internal energy of the fluid at the 
expense of kinetic energy of the turbulence. When these 
local losses due to dissipation are balanced by the local 
inputs of energy, then the turbulent structure in a given 
shear flow might be in a state of dynamical equilibrium. 
Since the turbulence consists of fairly large velocity 
fluctuations governed by nonlinear equations, one may expect 
a behavior like that exhibited by simple nonlinear systems 
Weeieelimit cycles. In a rapid energy transfer situations in 
which past events do not dominate the dynamics, one may 
expect that this limit cycle type of equilibrium is governed 
by the local length scales and time scales of the mean flow. 
In this turbulent motion the largest eddies may be as big as 
the width of the flow, which may then be an appropriate length 
scale in the analysis of the interaction of the turbulent and 
the mean flows. Large eddies lose a significant fraction of 
their kinetic energy within one "turnover time", but a 
negligible fraction is directly dissipated. This implies that 
the nonlinear mechanism that makes small eddies out of large 


ones is as dissipative as its characteristic time permits. 
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The nonlinear mechanism is ultimately dissipative because 
it creates smaller and smaller eddies until the eddy sizes 
become so small that viscous dissipation of their kinetic 
energy is dominant (eddy Reynolds number of order one). 
Since small-scale motions tend to have small time 
scales, one may assume that these motions are statistically 
independent of the relatively slow, large scale turbulence 
for a sufficiently large Reynolds number for the mean flow. 
Therefore, the small-scale motion should depend only on 
the rate at which the large-scale motion Supplies energy 
and on the kinematic viscosity. Thus it 1S assumed that 
this rate of energy supply should be equal to the rate of 
dissipation. If the dissipation rate, € , can be related 
to the length and velocity scales of the large scale tur- 
bulence, one can assume that the rate at which large eddies 
supply energy to small eddies is proportional to the 
reciprocal of the time scale of the large eddies and 
therefore viscous disSipation of energy can be estimated 
from the large-scale dynamics which do not involve viscosity. 
The amount of Kinetic energy per unit mass in the large- 
scale turbulence is proportional to eS and the eddy time 
scale is assumed to be proportional to ‘/ui?)“* where 1 
represents the size of the largest eddies or width of the 
flow. From this general idea, Garwood (1976,1977) derived 


= V2 
a total dissipation time scale, te = a + + aor 
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derivation see section I.D.1 and for more detail see 
Garwood (1977)) by including the planetary scale of 
dissipation. This is the integral time scale of the 
turbulence. Depending upon the time scale ( Cr ) of the 
surface forcing (variation in atmospheric conditions), an 


interaction between two time scales might be expected. 


B. SUBJECTIVE EXPECTATION AND ANALYTICAL SOLUTION TO 
INTERACTION BETWEEN GW AND © 


1. Constant Dissipation Time Scale 

Garwood (1976) developed a solution for forcing 
period much longer than the integral time scale of the 
turbulence. This solution is governed py the dissipation 
time scale (%). For these time scales one can assume a 
quasi-steady state so that (1-2) and (1-3) become diagnostic 
equations for gw> and ¢u®+¢+v?> . In this way, the 
entrainment buoyancy flux (Fig. 3) can be computed directly 
as a function of the surface fluxes of momentum and buoyancy. 
However, fluctuations in surface boundary conditions ae 
and bw'(o) of sufficiently short period require the con- 
Sideration of the unsteady terms. 

A very simple model which can be solved analytically 
Pee ides an indication of the importance of this unsteadiness 
and how it might be treated. In the simple model, only two 
large terms, shear production and dissipation, are balanced 


by the time rate change of the turbulent kinetic energy, 
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aE) , EX - Ww (iyeasinut ) (3-1) 
ot Ct hh 


where "a" is the amplitude of fluctuation and W= 2T/ Tr 
as its frequency. Assuming that the dissipation time 
scale is constant and defining a dimensionless time, Cats ee 4 
a dimensionless turbulent kinetic energy, p= BEEP ; 
and frequency, Wz=G@wWw2271%/G& 


(3-1) is written; 


a : 
= 1+ ASIN( We -2 
a + i +) (3-2) 


For cyclical steady state, the analytical solution to (3-2) 


is 





W=i+-( + ) C SiN Wat — We Cos Wet) (253) 
For the case a=o, YP =1 which is the non-cyclical 
steady state solution. For finite amplitude, a, the full 
solution (3-3) is needed. In terms of frequency, for the 

low frequency fluctuations, wy ¢<~-5 (i.e., long period 
compared with @ ), the approximation of quasi-steady 
state is satisfied. Depending upon the Rossby Number 

(Re = Us$h)) and the mixed layer stability parameter 

lag a a ), this low frequency fluctuation corresponds 


to a fiuctuation period of several hours or longer. From 
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the complete solution (3-3),the important aspect is the 
relative magnitude of the response associated with the 
fluctuating component (fig. 4). For very high frequencies, 
Ww >? 1, the amplitude of the response is negligible. In 
other words, the high frequencies are filtered out. This 
becomes a possible cause of error in using observed winds 
to drive a model having an integration time step smaller 
than or comparable with G& . If the quasi-steady state 
assumption is made to facilitate obtaining a solution, but 
without filtering or smoothing the surface boundary 
conditions, an incorrect high-frequency response will not 
only be present but it may bias the mean trend. 
2. Non-constant Dissipation Time Scale 

In this case the same simplified equation (3-1) is 
considered, but the dissipation time scale @ is no longer 
constrained to be constant. Using the convective dissi- 


pation time scale Gs wy of Garwood (1976), (3-1) becomes; 


=v 


3 


= =. 3/2 
= =- + 2 Cita sin wt) (3-4) 





To examine the significance of the unsteady term, equation 


(3-4) is first nondimensionalized by defining 


_ <> 
Wy uz ) 
Ta et , 
hy 
Wh 
and Uy, = ie 


By, 





then (3-4) is transformed to: 
ay 32 ; 
W +1 +a sin ( GT) 3-5 


The equation (3-5) is solved numerically for several 
Combinations at amplitude (a=0.1, 0.5 and 1), and angular 
frequency of forcing ( wh = 200m coe ,2c% |, Oe on, and 
@eeen). Initial adjustment was achieved at about vet 3 4 
There is also a tendency that the higher the forcing 
frequency, Wh , the shorter the initial adjustment period 
that is needed. At the early stage of the wind forcing, the 
change of turbulent kinetic energy 1s mostly due to the 
wind forcing itself. As the energy supply from the wind 
continues, and after an initial adjustment is achieved, the 
balance among the three terms 1S changed and becomes a 
function of the amplitude and frequency of the forcing. 
Different responses of the unsteady term 
to different frequencies of forcing are shown in figure 
6 a-b for the case of small (a=0.1) and large (a=1) 


amplitude fluctuations. For both amplitudes, there is a 


tendency for the mean value of the unsteady term, 4 , 


to vary with Preemie nee even though it is very small. fMThis 
1s seen more clearly in the case of small amplitude fluctu= 
Peetery (a=.1) because of the magnified scale for the small 
amplitude case. The mean value is a little above zero for 
high frequency and it approaches zero as the frequency is 


getting lower. This effect is also shown in figure 6 which 
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are plots of the case 1); SinWwt function of the wind 
forcing versus Wt , case 2); fluctuating component of 
dissipation normalized on the amplitude of the wind 
production ( (EP) /a ) versus wrt eieie =I, eral 
Case 3); same as case 2) except a=.l1. Comparing the 
menctuating component of dissipation, case 2) and 3), with 
the constant wind forcing case 1), for each of three fre- 
quencies, one can see that the higher the frequency, the 
smaller the fraction of wind production that is dissipated 
in the mean. In other words, more energy is stored in the 
mean to deepen the mixed layer. 


Figure 7 shows the peak-to-peak amplitude of the 


a ROLE? wh 
unsteady term us ot versus frequency iy for three 


emplitudes of fluctuation a=1,0.5and 0.1. For each case, the 
peak-to-peak amplitude is almost a constant close to 2,1 


and .2 respectively over the frequency range WA san ico 
am 


200 1 and it is negligibly small for the frequency 

range below .oO2TI » ome interesting transition 

cases occur over the frequency range Wh <= po2T to 211 
Ld 


in which the adjustment time scale and forcing time scale 

are comparable, and therefore interact. For relatively 

low-frequency surface forcing, Wh < .02zT (nondimensional- 
Li» 

ized period Ts = @M/te 2100), regardless of magnitude of 

oscillation amplitude, the turbulent kinetic energy budget 

is in a quasi-steady state, i.e., dissipation is nearly 


proportional to the wind shear production as indicated by 
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the negligible peak-to-peak amplitude of unsteady term, 


bh_3@> . Therefore, for this case, (3-4) reduces to 
up ot 
CEs uz (142 8in wt)? (Cua 


ferecelatively small amplitude oscillations (a <.0.5), 
regardless of frequency, it may be reasonable to assume 

a quaSi-steady state for many cases in which a high degree 
of accuracy is not required. But for relatively large 
amplitude oscillations, the unsteady term is more clearly 
frequency-dependent. For ad0.5 and in the transition 
frequency range below 28 =r lat eee MA yoo. Sumi. be 
a bad approximation to assume a quasi-steady state and use 
(3-6) in an approximate sense. On the other hand, for 
az2.5 and frequency who =. 2m , this unsteady term is 
no longer negligible. Therefore it is strongly suggested 
that the full form of equation (3-4) be used to avoid an 
incorrect high frequency response. 

Figure 8 shows the phase difference between the 
wind forcing and the dissipation as a function of the 
frequency of the forcing. For the frequency range 
there is almost no phase shift. As shown by the negligible 
amount of hysteresis in figure 9a, the dissipation is nearly 


proportional to the wind shear production after the initial 


adjustment period o< “ee <4 . In the other extreme of 
: : WwW : 
meen trequency forcing ( 4A 2? 207 ), there is almost a 90° 


phase shift between wind production and dissipation. As 
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shown in figure 9d, the adjustment period is longer 
relative to the forcing time scale, but dissipation 
converges to a constant value that is independent of the 
frequency and magnitude of the forcing oscillations. For 
the frequency range, .02T CWA £207 , in which the two 
time scales are comparable, the phase Eureenineeeae es 
nearly from zero at .0O2T to nearly 90° at 20T . Figure 
9b and 9c show that since the adjustment time scale and 
the forcing time scale are comparable and interact, this 
causes the ellipsoid to achieve a maximum width in this 
domain, 1<Te<10. An important result is revealed by 
figure 9b which shows that the amplitude of the response 
vere 1S nearly proportional to the amplitude of the 
forcing cycle, inspite of the hysteresis. Therefore, if 
mae@ce information 1s not crucial, it is accurate for this 


case to assume a quasi-steady state. 
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IV. ASYMPTOTIC PROPERTIES OF THER GARWOOD MODEL 


A. NUMERICAL RESULTS OF GARWOOD MODEL 

If quasi-steady state 1S not assumed, the full prognostic 
equations may still be solved numerically. The full 
equations of Garwood model, (2-5), (1-2), (1-3), (1-4) and (1-5) 
mere the jump conditions (1-6) and (1-7) are solved 
numerically for about a one-day period using hypothetical 
initial conditions, h=lm, ‘¢&>=.03 ¢m7/sec* , CW'2> =. 01 Cm7sec* , 
Do , Vso and three different cases of 4b (ab=5x16* , 
Exton? and 5x10? cm/sec? ). The size of 4b indicates the 
degree of stratification and is related to the lapse rate,T, 


as follows; 


_ SS Ge eer 
4b = > Lf bwW’‘(o) dt (4-1) 
ous oT Pn ne : 
where r= 28 = “355 > -25 Ss for the simple case of 


vertically homogeneous salinity. Assuming the surface 
buoyancy flux to be negligible, Pw (o=s0, Ab=5xI0% cm/sec? 


corresponds to Ss =1°/2z50m (i.e. f= 1075 sec ) ab=5xiotm/ee 


4 
corresponding to Sf=1%/as5m (i.e. P= 107+ secs? 2). ~+=42The 
constant employed here is mz@€ . To simplify the problem, 


the momentum flux U’w’(o) =I Cm?/’sec?, v’W'so0 , the constant 
frictional velocity Uy= 1 ¢m/sec ,Qs0, Pw'so and f=i0+ sec’ 
are also assumed. Using these conditions, three different 


curves of each parameter are plotted for each lapse rate. 
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Figure 10 shows that the stronger the stratification, the 
larger the amount of turbulent kinetic energy in the mixed 
layer. It starts from ¢ E)=.03c¢mecand very rapidly 
increases during a short period immediately following an 
impulsively applied surface forces. £E> shows a maximum 
at t=3~5 hrs., a minimum at t=17~18 hrs. and a second 
Maximum at t=24~25 hrs. The figure also shows that the 
weaker the stratification, the earlier the maximum is 
reached. The vertical component of turbulent kinetic 
energy shows a Similar pattern except for the case of 
weak stratification (= =1*c/as0m) in which it decreases 
monotonically (fig. 10b). Figure 11 shows <w°>®) versus 
mame for <r =1°C/2.5m »I%/2em and \°C/agom +. Starting with 
iw y@= 0.33, it reaches 0.15, 0.14 and 0.13 respectively 
after a one-day period. These sudden reductions in <W’*>£é> 
occur at a very early stage and then stay nearly constant. 
Bamounier words, vertical components of turbulent kinetic 
energy relative to total kinetic energy adjust rapidly in 
the early stages, and, the stronger the stratification, the 
larger will be the vertical component relative to total 
kinetic energy. 

Figure 12 shows changes in Ab with time. Each ab 
increases monotonically, although more rapidly in the case 
of stronger lapse rate. Larger values of Ab contribute to 


more turbulent damping due to mixing at the interface. 
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Figure 13 shows the change in current with time, and 
indicates that the stronger the stability, the larger will 
be its magnitude. With the imposed east wind, the currents 
start flowing westward at the initial instant. The west 
component reaches a maximum value of 17 cm/sec, with the 
mMemthward component current still growing. As the west 
component decreases to zero, the north component current 
reaches a maximum of 25 cm/sec. This slowly turning current 
Seerection is caused by Coriolis force which deflects the 
current direction to the right in the Northern Hemisphere. 
After one inertal period, both GO and V decrease to zero. 

The maximum of U*+¥* occurs at about t=4hrs and the 

Minimum is at about t=l17hrs. Also the figure 13 shows 

that the magnitude of each successive maximum value decreases 
with time. 

Figure 14 shows the monotonically increasing mixed layer 


depth with time. All cases start from lm and after 10sec 


the mixed layer depth is 45m for 23 = 1°°/250m, Zam 1 Or 
T 
<= = oC 725m anieusalim Tow Fey = 1S fe oie Theses ine) 


cate that the larger the stability, the lesser the deepening 
mo iS. 
: a d Wray dh dh 
Figure 15 a-c show on /<w 2 ; Se [Us (or ais 
because U,=lcmMeris assumed) and di fe >¥? for each lapse 
rate. The middle curve of each figure indicates that 


decreases more rapidly in large stability than in weak 


e e dh 
Seaoklity. F a ape Se zyVvz_sCS, d Y. 
a “L1G or each case on) ah / Uy an dh, /< 8 > 
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Pee notebarallel prior to about one inertial period, li.e., 
BH ce Cw and oN ce Us would not produce the same deepening 
even though Bw'(e)so is assumed. The figures imply that in 
general the deepening is not proportional to Lw'>¥? , 467 


Sete, OUL rather depends upon the full form of entrainment 
dh _ ce<Wv2 
a ~ > 4+hSb 
SO that gh =¢wre? during the early stages (Figure 15d). 


eamation An exception occurs if hab<=o 
Figure 16 shows the inverse of the bulk Richardson 
number changing with time. For t <20 mins, R; is large 
in the case of weak initial stability because of the small 
Ob. After 20 mins, the larger values of Ri? occurs for 
the cases with larger initial stability. This is due to the 
strong mean shear developed at this stage. At one inertial 
period, Rim goes to zero because of the stagnation in 
current at this moment. Halpern sii Suggested that if 


local Richardson number (Rn = 22 / 4 ) is less than 


aa 
Sea ourthn, instability may occur, and if Ry =1, marginal 
Stability exists. He also suggested in his experimental 
study of Norther Pacific Ocean that QR.) =z is a critical 
Semadition for the occurence of shear instability. Therefore, 
for the case of £1 =1°c/a5om, there is an instability at 
meowe t=lLOOsec. If Rit 1s considered to be an average of 
Ry over the mixed layer depth and therefore Ry > Ran ; 
there may also be an instability at t=4hrs for the case of 


oT = I°C/a.ism - This instability leads to a maximum turbulent 


kinetic energy at this time interval. 
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Therefore, how the turbulence develops with time, as 
indicated by equation (1-2) and (1-3), can be summarized. 
fier total turbulent kinetic energy (fig. 10) 2—&> jumps 
due to impulsively applying surface stress to an initially 
resting fluid, whereas the energy produced by any mechanism 
is dissipated at a rate proportional to LE > ee Woksvaue Hans, 
shear production is increased to a maximum value even 
though the damping also increases. The maximum turbulent 
kiretic energy exists at about t=4hrs, and then decreases 
toward a minimum mainly due to decreasing shear production 
at about one inertial period. 

The more detailed, term-by-term study is left for the 


mext section. 


-B. ASYMPTOTIC REGIMES IN MIXED LAYER DEEPENING 
Combining equations (1-2) and (1-3) of the Garwood 
model, the total turbulent kinetic energy budget equation 


1s obtained. 
dc hee») { = 2) - 2m Us + jazt*dh + Wh Bw(o) ~hab = ~2 EY? 4 £hce)) 


Using hab=tn'h? , one can write the above equation in the 


following order; 


tyth2sh ~|az}2 2 ie ne = 2mu; -2 (ss + $hdE>) (4-2) 


(A) (B) ec) (D? 


The corresponding equation from de Szoeke and Rhines (1976) 


may be obtained by dividing equation (1-30) by a 
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4 
2 va 24dh — s 
[ ¢8*h _ rie (\- cos $t) + G Us ae = Z2M,Us (4-3) 
(A) (8) CC) Uy? 


Equation (4-3) is originally the Niiler model except for 
the term (C). The term(A)in both (4-2) and (4-3) is an 
entrainment damping term which is a upward buoyancy flux 
Paeune Dase of this mixed layer. The term (B) in both 
equations is the shear production term. The term (C) 
in (4-3) is the energy needed to spin up the level of 
turbulence, Cus , in the increment of mixed layer dh in 
memerdt, while the term (C) in (4-2) is the nonstationarity 
term which shows the incremental change of the turbulent 
energy budget d(h¢&>) in time dt. The term (D) in both 
equations is the wind shear production minus dissipation. 
Before analyzing the numerical solutions it might be useful 
to examine the characteristics of both equations and their 
differences. 

men bouy models Contain nonlinear equations, i.e., the 
Garwood model is a set of six nonlinear differential 
equations which can be solved implicitly for each parameter, 
while the de Szoeke and Rhines model has one nonlinear 
differential equation. 

ii) Both models have numerical problems when the zero- 
mmebtoa?! condition (h(t=0)=0) are used, although such 


eem@artions are not likely to occur in nature. 
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iii) The nonstationarity term in Garwood, - 
n 46S 4By dh is not exactly comparable to Ce, Up th. 
of de Szoeke and Rhines. Even with the assumption of no 
surface buoyancy flux, <EYuy does not approach a constant 
ino 1 = \%/2.5m, ST = [°c/25 m and 4T = 1° /250m Lipo 
one inertial period. After this period, Gu, = constant 
holds approximately, although it is still not exact except 
Mer the case of weak stratification (figure 10a). Even 
when the condition <EyYu, “= constant holds, there is a 
difference of ioe between the two models. 

fe) erm (D) an the Niiller (1975) model is 


parameterized in the manner; 


aw’(2.E _ ere All eee Oe = 
w(£+8)I | dene z (édz M. Ux 





where "d" is the depth of the perturbation energy prod- 
uction zone near the surface. This implies that only the 
near surface production is dissipated. In other words, 
dissipation is a fraction of production term near the 
surface. None of the shear production energy at the base 
of mixed layer and the surface buoyancy production energy 
(this term is not considered for this study) goes to 
dissipation. Because of this parameterization, the shear 
production at the base of mixed layer has only an interaction 
with the entrainment damping term. On the other hand, in 
the Garwood model, dissipation is parameterized in a way 
that all the production including the entrainment shear 


production and the surface buoyancy flux is dissipated. 
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Term (D) in the Garwood model is the turbulent production 
myer the wnole depth minus the total dissipation, i.e., 
the viscous dissipation as well as the planetary scale 
dissipation. This second dissipation term is not included 
in the Niiler model. If the planetary scale dissipation 
LE? 
Us 
term (D) in the Garwood model approaches term (D) in the 


is neglected and = a constant is considered, then 
de Szoeke and Rhines equation of the Niiler model within 
a constant factor. 

With these characteristics and differences in mind, two 
cases of the Garwood model and one case of the de Szoeke and 
Rhines model are solved numerically using the same conditions. 

Case 1) The equations (2-5), (1-2), (1-3), (1-4)and (1-5) 
with the jump conditions are solved numerically using the 
Same constants and initial conditions as in the case of 
QT = |°c/25m (fsio*sec™* ) in Section IV.1 (i.e., 
Parameters and boundary conditions: M=€, Un=1n/sec , $= 1o~* see! 
Uw! (0) =1 em3/sece2 , yw’ o)20, YW'o)z=0 and initial conditions: 
h=loocm, <E>=2.03¢m*/sec2 » <W>=.0] cm2/sec? » 4625x107 cm/sec?) 

@ase 2) is similar to case 1) except that there is no 
phanetary scale dissipation term in (1-2) or (1-3). 

Case 3) Equation (4-3) of de Szoeke and Rhines is solved 
numerically using G=m,=! » Ugsl cm/sec » £= io 4 sect 
and M=N?s |o+ sec 

Fach solution is displayed from t=50sec because the 
iaeral values of terms (C) and (D) are large. The three 


cases of term (A) are plotted in Figure 17. Maximum 
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entrainment damping occurs at t=4~5 hrs and a second 
maximum at t=25 hrs for each case. The entrainment damping 
of the de Szoeke and Rhines model is twice as large as that 
of the Garwood model at the time of the first maximum and 
remains larger than that of the Garwood model throughout. 
@ese 2) of the Garwood model and Case 3) of de-Szoeke and 
Rhines are very close at time equal to one inertial period. 
The entrainment damping of Case 2) (i.e., without planet- 
ary scale dissipation) is larger than that of Case 1) of 
Garwood model after t=15 naoand the difference increases 
with time. For the time interval less than 15 mins, the 
planetary scale dissipation term does not make any con- 
tribution to the entrainment damping. The large damping of 
term (A) in the de Szoeke and Rhines equation compared to 
the Garwood model occurs inspite of the large shear prod- 
uction in the Garwood Model (Figure 18). This may be due 
Gempne Tact that none of the shear production of the 

de Szoeke and Rhines equation dissipates. Instead, all 

of this energy goes to deepen the mixed layer, and 
therefore the entrainment damping grows larger to offset 
iemratser entraining energy. The reason for larger 
damping in case 2) than in case 1) is because there is less 
dissipation in case 2) (without planetary scale dissipation) 
and thus there is more entraining energy available, against 
which entrainment damping grows. Since this planetary scale 
dissipation increases with time, the difference in the 


entrainment damping betiveen the two cases also increases 
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with time. Since the entraining due to entrainment shear 
production decreases to zero (Figure 18) at one inertial 
period in all cases, the net entraining energy in case 2) 
and case 3) is the same at this time stage and therefore 
the entrainment damping in case 2) and case 3) is almost 
the same. Whereas, case 1) has less entraining energy 
because of planetary scale dissipation. This leads to less 
entrainment damping at this same time stage. 

Figure 18 shows term (B) for the three cases. The 
figure shows that the entrainment shear production of the 
Garwood model is not affected by the planetary scale 
dissipation term and the term (B) of the Garwood model 
is larger than that of the de Szoeke and Rhines model at 
least for one-day period. It also shows that the time of 
maximum shear production agrees with the time of maximum 
entrainment damping as discussed in term (A). The reason 
for the smaller entrainment shear production of the de Szoeke 
and Rhines model is because there is no dissipation 
attributable to the entrainment shear production, and this 
leads to larger deepening in the de Szoeke and Rhines model 
than in the Garwood model (Figure 19). Therefore, the term 
(B) of the de Szoeke and Rhines model, which is divided by né 
Shows a smaller value of extrainment shear production. 

Figure 20 shows the value of the (C) terms. Again, the 
planetary scale dissipation term of the Garwood model does 


feet tect this term. This term in the Garwood model is 


DS 
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larger than that of the de Szoeke and Rhines, although 

the ference between the two models diminishes in time. 
The reason for term (C) (Figure 20) being larger in the case 
of the Garwood model, and possible for the decrease in the 


difference with time is because of the omission of the term, 


hs » in the de Szoeke and Rhines model. . In addition, 
this is also because LE>/ Ux tends to be large at first and 


only becomes nearly constant after about an inertial period. 
iteetasc and important (BD) terms are plotted in figure 
Pie The term (D) of de Szoeke and Rhines is equal to 2.0 
and thus appears as a Straight line. On the other hand, 
the term (D) for case 1) and case 2) of the Garwood model 
varies with time. After several minutes, the planetary 
scale dissipation term starts to act and because of this 
getion, the term (D) of case 1) is always less than that of 
case 2) with larger differences at the maximum and minimum 
values. After one inertial period,case 2) and case 3) are 
close together, because there is no entrainment shear, 
therefore no dissipation of entrainment shear production 
(this is always the case for de Szoeke and Rhines) and no 
dissipation from planetary rotation in both cases. The 
larger differences of the term (D) between case 1) and case 
2) at the time of minimum is because of the larger value 
of <E> due to large entrainment shear production which 
contributes to larger planetary scale dissipation, 


The term (D) of the Garwood model shows a negetive value 
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between t=l.3 hrs and t=10.1 hrs. The reason for this 
negative value during this time interval is because the 
fermm (D) is not a net production term, i.e., this 
term does not include the shear production at the base of 
the mixed layer but does include the dissipation from that 
entrainment shear production. Therefore, this term shows 
a Minimum during this time interval whereas term (B) 
(Figure 18) shows its maximum. 

To examine the asymptotic regimes, the solutions for 
case 1), case 2) and case 3) are plotted for separate 
may, (8B), (C), and (D) terms and also for B/A, C/A and 
D/A following de Szoeke and Rhines (see also Figure 2). 
However, Since the results for case 1) and case 2) are al- 
most the same, only case 1) will be considered. From the 
moectal instant to about t=5 mins a balance (C)~(D) holds 


for the Garwood model (Figure 22). 


se ~ 2muU2 -2 (EM 4 Sher) (4-4) 


This is similar to the balance in the de Szoeke and Rhines 
model except that the duration for their model is less than 


100 sec. In the de Szoeke and Rhines model 


Co ug go me 1M, UR Or dh oc Uy Gres) 


2 





Mas (C)~—(D) balance for both models implies that the 
initial deepening proceeds due to turbulent motion in the 
mixed layer. 

After this period, approximately from t=5 mins to 
t=18 mins, there is no simple balance in the Garwood 
solution. Therefore, the full form of the turbulent 
kinetic energy equation is required during this period. 
However, in the de Szoeke and Rhines model, a balance 
menween (A) and (D) holds rather clearly from the early 
stage to about 1.2 hrs. For the Garwood model, this is 
the period that the value of nonstationarity is decreasing 
and entrainment shear production is increasing, although 
both the entrainment damping and the production minus 
dissipation are increasing to offset the value of non- 
stationarity and shear production. From t=18 mins to 
about t=l12hrs, a balance (A)-(B) holds for the Garwood 


model, 
1 2ap2zdh 2 dh 
meg, ~ (4c Fe 
or 
2 
= Leh HK (4-6) 
hab 
This is the regime in which the bulk gradient Richardson 
numbers is nearly constant (R;z! ). For this time domain, 
the value of nonstationarity is decreasing and the strong 
shear production is balanced by the entrainment damping. 


In the de Szoeke and Rhines model, this balance (A)~(B) 
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elds from t=1.2 hrs to 8.8 hrs, which is somewhat shorter 
period compared to the Garwood model. Finally following 


t=12 hrs, a balance (A)~(D) holds for the Garwood model 


Suth? 4h Pehl: Erde ence (4-7) 


Likewise in the de Szoeke and Rhines model, the same 
balance (A)~(D) holds although this begins earlier at 


m—o.o hrs. 
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V. SUMMARY 


The one-dimensional model of the ocean mixed layer of 
Garwood (1976, 1977) has been further generalized. The 
entrainment buoyancy flux equation yal welalsutss anroyel’S IL abysyfelisrezb vtec 
from the turbulent kinetic energy budget at the base of the 
mixed layer. In earlier models, the unsteady term 
has been considered to be negligible in the derivation of 
the entrainment buoyancy flux equation. From the reduced 
form of the turbulent kinetic energy equation in which a 
balance holds among the shear production, the dissipation 
and the unsteady terms, Garwood (1976) nondimensionalized 
andsolved this equation linearly using a constant dissipa- 
tion time scale. He suggested that if the quasi-steady 
state assumption is made to facilitate solution, without 
filtering the surface boundary conditions, an incorrect 
high frequency response will not only be present but will 
bias the mean trend. Following these suggestions, the 
unsteady term has been included in the parameterization of 
entrainment buoyancy flux equation. Even though the entrain- 
ment buoyancy flux equation of the Garwood model is sufficient 
in most cases, this new generalized equation is necessary 
to study the asymptotic regimes. 

To examine the importance of the unsteady term in the 
kinetic energy equation in more detail, the same reduced form 


as the Garwood (1976) equation is nondimensionalized, and 
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solved numerically using a nonconstant dissipation time 
scale, %=h/<€>* . Solutions for variable amplitudes of 
fluctuation and forcing frequencies were obtained and 
their results are as follows. After an initial adjustment 
period od wee 4 

i) For any amplitude of fluctuations, the quasi-steady 
state assumption may be good, ifWh/u, €0.02T. 

mijewein the case where the amplitude of the fluctuation 
/mememaller than a half of the mean amplitude or Wh/UW,< .2T 
(i.e., Te> 10) the assumption of quasi-steady state may be 
possible. 

iii) For the frequency Wh/Use>.2T1 , or the amplitude 
perme larger than one-half of the mean amplitude, the 
unsteady term should not be neglected. 

Following the lead of de Szoeke and Rhines (1976) and 
fember (1975), the relative contributions of entrainment 
damping, entrainment shear production, nonstationarity, 
wind shear production and dissipation were examined. The 
Specific differences between two models are: 

i) The Garwood model is computed implicitly for each 


term whereas , the de Szoeke and Rhines' model is solved 


Sapereitly. 
ii) Even with assumption.of <B/us = a constant, an 
additional term he’ is in the Garwood model for the non- 


Seaceonarity term (C). 
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iii) Dissipation is parameterized in a different 
manner, that is, Niiler (1975) parameterized the 
dissipation term as a fraction of near-surface production 
term, whereas Garwood (1976, 1977) parameterized the 
dissipation term in relating to net production. The 
numerical solutions of both models indicate that larger 
entrainment damping occurs in the de Szoeke and Rhines 
case inspite of smaller shear production compared to Garwood. 
Also, the solutions indicate that a larger nonstationarity 
term occurs in the Garwood model, and that a constant 
production minus dissipation term occurs in the de Szoeke 
and Rhines model compared to large variation in the 
Garwood model including the negative value of this term. 

Following the method of de Szoeke and Rhines, each 
asymptotic regime is classified for both of the models. 

imeerom the initial instant to about 5 mins. a balance 
between the nonstationarity and the production minus 
dissipation holds in the Garwood model, i.e., 


SS Ee 2k —2 (ES 74 $4 <¢E)) 


(C) (D) 
This same balance holds in de Szoeke and Rhines with 
duration from the initial instant to about less than 100 sec. 
ii) From t=5 mins to t=18 mins, a balance holds for the 
mee torm ot the turbulent kinetic energy equation, whereas 


in the de Szoeke and Rhines model, a balance between the 
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entrainment damping and the entrainment shear production 


holds; 
tenth’ fh ES jazi* da 
(A) (B) 
or 212 
_. [4c 
Ri = Ab 


This balance (A)~(B) holds also in the de Szoeke and Rhines 


Pade rrom t=l_.2 hrs to t=8.8 hrs. 


io following t=12 hrs, a balance between the entrainment 


damping and the production minus dissipation holds; 


212 = 2/ = 
5N7h aA ~ 2mu -2 (ce>*+ $he>) 
(A) (B) 


fatemess in the de ozoeke and Rhines model, the same balance 
(Ay~(D) holds, although this balance occurs a little 


eewanarer, from t=8.8 hrs. 
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Figure 1. Idealized model for ocean mixed layer (----). Mixed layer depth 
is (h};(§) is the thickness of the interface or entrainment zone. 


(Garwood (]976)) 
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Figure 2. Numerical solution of the complete mixed 
layer model equation of ce Szoeke and 
Rhines showing various esymptotic regimes. 
Solid curve: the ratio D/A of terms in the 
equation; dashed curve: B/A; dotted curve; 
C/A. Parameters u,=l1,ms,=l,c,=1 (de Szoeke 


and Rhines(]976)) 
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figure 3. General solution to entrainment and turbulent kinetic energy 


equations. (Garwood (]976) ) 
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Figure 7. Peak to peak amplitude of unsteady term 


vs. nondimensionalized frequency. 
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Figure J]7. Comparison of the entrainment damping 


term in three cases 
Case 1) Garwood--full equations, 


Case 2) Garwood--without planetary scale dissipation. 
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Figure 24 
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